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Prediction and Validation of Mars Pathfinder Hypersonic 

Aerodynamic Data Base 
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Iv. James Weilmuenster, Robert A. Mitclieltree, 

Walter C. Engelund, Richard W. Powell 
NASA Langley Research Center 
Hampton, VA 23681-0001 

Postflight analysis of the Mars Pathfinder hypersonic, continuum aerodynamic data 
base is presented. Measured data include accelerations along the body axis and axis 
normal directions. Comparisons of preflight simulation and measurements show good 
agreement. The prediction of two static instabilities associated with movement of the 
sonic line from the shoulder to the nose and back was confirmed by measured normal 
accelerations. Reconstruction of atmospheric density during entry has an uncertainty 
directly proportional to the uncertainty in the predicted axial coefficient. The sensitivity 
of the moment coefficient to freestream density, kinetic models and center- of-gravity 
location are examined to provide additional consistency checks of the simulation with 
flight data. The atmospheric density as derived from axial coefficient and measured 
axial accelerations falls within the range required for sonic line shift and static stability 
transition as independently determined from normal accelerations. 


Nomenclature Introduction 


Aa axial acceleration, earth g 

An normal acceleration, earth g 

Ca axial force coef. 

Cd drag coef. 

C„, pitching moment coef. 

C m ,a pitching moment coef. 

derivative, per radian 
Cn normal force coef. 

D probe diameter, m 

h altitude, km 

M Mach number 

m entry mass, kg 

p pressure, N/m 2 

T© freestream temperature, K 

S reference area, m 2 

Voo freestream velocity, m/s 

x, y,z Cartesian coordinates, m 

a angle-of-attaek, degrees 

ax total angle-of-attack 

(a 2 + / 3 2 ) 1 / 2 , degrees 
3 angle of yaw, degrees 

7 effective ratio of specific heats 

Poo freestream density, kg/m 3 
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T HE hypersonic aerodynamic data base used in 
the Mars Pathfinder entry, descent, and landing 
(EDL) was predominantly generated using computa- 
tional fluid dynamics. While ground-based and flight 
data from the Viking project were used to check por- 
tions of the aerodynamic data for Pathfinder, differ- 
ences in entry velocity and angle-of-attack between 
Pathfinder and Viking required a more critical eval- 
uation of high temperature effects on hypersonic aero- 
dynamics. A systematic evaluation of Pathfinder aero- 
dynamics as a function of velocity, altitude, and angle- 
of-attack in a continuum regime revealed two intervals 
during a nominal entry in which static aerodynamic in- 
stabilities would induce vehicle “wobble”. The static 
aerodynamic instability was associated with changing 
sonic line location and associated pressure distribu- 
tions as a function of chemical kinetics and flow energy. 
These predictions and explanations were published in 
journal articles in 1996 1,2 and the “wobble” was con- 
firmed by Pathfinder flight data (accelerometers) on 
July 4, 1997. The physical phenomena leading to 
the static aerodynamic instability are reviewed, the 
measured accelerations confirming the instability are 
discussed, and additional post-flight analyses are pre- 
sented. The fidelity of comparison between pre-flight 
prediction and flight data significantly reduces uncer- 
tainty in the atmospheric reconstruction process and 
provides enhanced confidence in the application of 
CFD to future planetary missions with reduced design 
margin. 
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Fig. 1 Diagram of the Mars Pathfinder probe. 

Preflight Analysis 

The Mars Pathfinder configuration is shown in Fig. 
1. The forebody is a 70° half angle cone with base 
diameter 2.65 m, nose radius 0.6638m and shoulder 
radius 0.06625 m. The on-axis reference point for aero- 
dynamic moment coefficients is located 0.662 m behind 
the nose. The actual center-of-gravity (CG) location 
is (xcG,ycG,z C G) = (0.000107, -0.000435, 0.71541) 
where s is distance behind the outer mold line of the 
nose and x and y are orthogonal to the body axis. 

Predictions of a static aerodynamic instability for 
the Mars Pathfinder vehicle associated with shifting of 
the sonic line from the shoulder to the nose and back 
have been described in the literature . 1,2 Specific con- 
clusions of the first citation relating to the occurrence 
of the static aerodynamic instability on the 70 degree, 
spherically capped Pathfinder configuration are sum- 
marized here. 

(1) As the Mars Pathfinder probe descends through 
the Mars atmosphere the minimum value of the post- 
shock effective 7 first decreases from frozen gas chem- 
istry values to equilibrium values (1.094) correspond- 
ing to a velocity of 4.86 km/s. As the probe continues 
to decelerate through a near equilibrium post-shock 
gas chemistry regime, the value of 7 increases, until 
reaching its perfect gas value at parachute deployment 
(0.42 km/s). 

(2) For the spherically blunted, 70 degree half-angle 
cone at small angles of attack the sonic line location 
shifts from the shoulder to the nose cap and back again 
on the leeside symmetry plane because of the change 
in 7 . This transition in sonic line location is shown in 
Fig. 2 for a 2 degree angle-of-attack. The subsonic flow 
region, as defined by the frozen sound speed, appears 
as the light gray region. The darker region represents 
zones of supersonic flow behind the bow shock. At a 
Mach number 22.3 on the preflight Pathfinder simula- 
tion, the sonic line on the leeside dips deep into the 
boundary layer over the spherical nose cap. Later in 
the trajectory, at Mach 16.0, a subsonic bubble forms 



Fig. 2 Sonic line location over the Mars Pathfinder 
at 2 deg. angle-of-attack and Mach numbers equal 
to 22.3 (left), 16.0 (center), and 9.4 (right) showing 
the effect of gas chemistry. 



x, m 

Fig. 3 Pressure distributions in the plane of sym- 
metry over the Mars Pathfinder at 2 deg. angle-of- 
attack and Mach numbers equal to 22.3, 16.0, and 
9.4 showing the effect of gas chemistry. 

immediately behind the bow shock just downstream of 
an inflection point in the bow shock. The sonic line 
from the nose skirts above the edge of the boundary 
layer. Still later, at Mach 9.4, the bubble expands to 
merge with the subsonic region at the boundary layer 
edge. The inflection point on the bow shock disap- 
pears. 

(3) Pressure distributions on the cone fustrum ap- 
proaching the shoulder tend to be very flat when the 
sonic line sits forward over the spherical nose, as shown 
in Fig. 3 for Mach 22.3. Effects of the expansion 
over the shoulder can only be communicated upstream 
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through the subsonic portion of the boundary layer. In 
contrast, pressure distributions on the cone fustrum 
approaching the shoulder tend to be more rounded 
when the sonic line sits on the shoulder, exhibiting 
a more pronounced influence of the expansion on the 
upstream flow as shown in Fig. 3 for Mach 9.4. 

(4) In general, windside pressures exceed leeside 
pressures on the cone fustrum, producing a stabi- 
lizing moment that pitches the probe back to zero 
angle-of-attack. However, the behavior of the pres- 
sure distribution in the vicinity of the shoulder sig- 
nificantly influences the pitching moment coefficient 
because of the larger moment arms and surface area 
at the edge of the probe as compared to the inboard 
nose and fustrum regions. For the Mars Pathfinder 
probe at 2 deg. angle-of-attack, the flat, leeside pres- 
sures approaching the shoulder (when the sonic line 
sits over the nose) can exceed the rounded windside 
pressures approaching the shoulder (when the sonic 
line sits over the shoulder). The net effect of this 
crossover distribution near the shoulder tends to pitch 
the probe to higher angles of attack (static instability). 
The overall balance (crossover-point location) is sen- 
sitive to freestream conditions and corresponding gas 
chemistry within the shock layer. The combination of 
small angle-of-attack, large cone angle, windside sonic 
line over the shoulder, and leeside sonic line over the 
nose with a secondary, subsonic zone above the leeside 
fustrum yield conditions which favor (but do not nec- 
essarily guarantee) a destabilizing pitching moment. 

(5) Representative parts of the aerodynamic data 
base as a function of velocity during the hypersonic, 
continuum portion of the trajectory are shown in Fig. 

4. Conditions for a positive, destabilizing moment 
coefficient derivative C m ^ a occur twice in the Mars 
Pathfinder mission. (C m is zero when ax equals zero; 
Cm, a is positive where the C m curve for a = 2 deg. 
crosses the C m = 0 axis in Fig. 4 b.) The first occur- 
rence (7.5 > Vqo > 6.5 km/s, 51 > h > 37 km, vicinity 
of peak heating for this trajectory) results from the 
transition in the sonic line location as a function of 
gas chemistry changing from nonequilibrium to equi- 
librium. The second occurrence (4.0 > V/g > 3.1 
km/s, 25 > h > 22 km) results from the transition in 
the sonic line location as a function of decreasing flow 
enthalpy in a near equilibrium gas chemistry regime. 

Measured Accelerations 

Measured accelerations during Pathfinder’s entry, 
descent, and landing (EDL) validate these earlier con- 
clusions. The raw, three-axis accelerometer data for 
the hypersonic portion of EDL are presented in Fig. 

5. In this figure, the z-axis is in the axial direction. 
The x- and y- axes are fixed on the body as it spins 
(approximately 2 rpm), pitches and precesses during 
descent. Peak dynamic pressure occurs approximately 
78 seconds after entry interface. In the vicinity of peak 



2000 4000 6000 8000 

V , m/s 


a) Axial coefficient 



2000 4000 6000 8000 

V , m/s 


b) Pitching moment coefficient 

Fig. 4 Predicted aerodynamic coefficients as 
function of velocity and angle-of-attack for early 
Mars Pathfinder trajectory with ballistic coefficient 
equal to 45. 

dynamic pressure, the magnitude of the off-axis accel- 
erations abruptly decrease, corresponding to a brief 
period where CFD simulations show the sonic line 
fixed to the spherical nose at small ax- 

The trend is more clearly observed in Fig. 6 in which 
a total normal acceleration (An = (.4^ + Ay) 1 /' 2 ) is 
plotted as a function of time. In this figure, a rather 
abrupt increase in normal acceleration is observed on 
either side of the relative low level (An ~ 0.0 14g) 
near peak dynamic pressure to a relatively high level 
(An ~ 0.04g). This change closely corresponds to the 
predicted times when the sonic line shifts from the 
shoulder to the nose (f» 71s) and then back again to 
the shoulder (« 82s). 
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Fig. 5 Science accelerometer measurements 
recorded on Mars Pathfinder during 2.5 minutes 
following entry interface with the Mars atmo- 
sphere. 



Fig. 6 Total off-axis acceleration recorded dur- 
ing 100 second interval surrounding peak dynamic 
pressure. 

In general, the correspondence of measured accelera- 
tions and predicted aerodynamic coefficients is compli- 
cated by the unknown freestream density. Velocity can 
be derived from integrated accelerations but density is 
not directly measured. This complicating factor can 
be removed by examining the ratio of normal to axial 
forces and normal to axial aerodynamic coefficients. 

^ F N t 4 

Cn _ i/- 2 Poa v'^s _ mAjsi _ An 

Ca ~ 1/2 p F J V *S ~ ™ Aa ~ Aa 

The predicted values of normal to axial force coef- 
ficient using the initial state vector determined four 
hours prior to entry interface in the six-degree-of- 


freedom (6-DOF) Program to Optimize Simulated 
Trajectories (POST) 3 code and aerodynamic coeffi- 
cient tables derived from CFD simulations 1 are pre- 
sented in Fig. 7a. The corresponding measured values 
of normal to axial accelerations are presented in Fig. 
7b. The predicted and measured results show very 
good qualitative and quantitative agreement. For ex- 
ample, the effects of the predicted and measured static 
aerodynamic instability between 82 and 90 seconds are 
evident in both the magnitude and duration of the 
pulse. The time averaged aerodynamic coefficient ratio 
in the vicinty of peak dynamic pressure between 71 and 
82 seconds is also in good agreement with magnitude 
and duration of the measured ratio of accelerations 
when the sonic line sits over the spherical nose. The 
predictions show a much larger oscillation about the 
mean in this domain than was measured. An overly 
conservative specification for C m<q in this region in the 
simulation is the suspected cause. The magnitude and 
frequency of oscillations beyond 90 seconds are in fair 
agreement but the mean value of measured accelera- 
tion ratio (.0016) remains at a constant level while the 
mean value of the simulated acceleration ratio starts 
lower (.0008) and gradually rises to .0024. Reasons for 
these discrepancies are not yet understood. 

The total angle-of- attack ar predicted in the 6-DOF 
POST analysis is presented in Fig. 8a. The corre- 
sponding values in flight (Fig. 8b) are derived from 
aerodynamic tables of Cn/Ca as a function of ax 
using the measured ratio of normal to axial acceler- 
ations. Here again, both qualitative and quantitative 
agreement are generally good and serve to validate the 
simulation methodology for aerodynamics. 

Atmospheric Reconstruction 

Reconstruction of the Mars atmosphere based on 
measured accelerations on Mars Pathfinder has been 
described by Spencer et. al. 4 Velocities are derived 
from measured accelerations. Densities are derived 
from measured axial accelerations, the derived veloc- 
ities, and the predicted aerodynamic data base for 
Ca as a function of velocity, Mach number, and total 
angle-of-attack. The point-by-point transformation is 
given by 

Errors in the predicted aerodynamic data base are lin- 
early related to errors in the derived densities. The 
reconstructed density is plotted as a function of veloc- 
ity in Fig. 9. The off-axis accelerations are plotted in 
this same figure to show the correspondence of simula- 
tion coordinates (velocity, density) with the onset and 
decay of the aerodynamic instabilities at small angles 
of attack. 
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t, S 

a) Predicted ar from 6-DOF POST 



t, s 

b) Measured normal to axial acceleration ratio 

Fig. 7 Ratio of normal acceleration to axial accel- 
eration recorded during 100 second interval sur- 
rounding peak dynamic pressure. The ratio is 
equivalent to the value of Cn/Ca for Pathfinder at 
its total angle-of-attack aj. 

Uncertainty Analysis 

The uncertainty in center-of-gravity (CG) location 
is ±5mm in the axial and ±1 mm in the radial direc- 
tions. Given Ca ~ 1.7 over the hypersonic, continuum 
portion of the trajectory, the uncertainty in C m asso- 
ciated with the uncertainty of the CG location in the 
pitch plane is A C m = C/jAr/T* = ±0.00064. The un- 
certainty in C m associated with an axial translation of 
the CG is two orders of magnitude smaller for small 
angles of attack. Note that there is no uncertainty in 
the aerodynamic reference point; we simply treat CG 
location uncertainty as equivalent to a corresponding 



t, s 

b) Derived ut from recorded accelerations 

Fig. 8 Comparison of total angle-of-attack cut pre- 
dicted by 6-DOF POST using initial state vector 
determined 4 hours prior to entry interface with 
values derived from measured accelerations. 

change in the reference location. 

All of the solutions generated for the next section 
and most of the solutions used to generate the aerody- 
namic data base were obtained on a surface grid of 30 
by 60 cells (half body to symmetry plane) with 64 cell 
normal to the body. A grid refinement to 40 by 80 cells 
for a case with \ ). x . = 3515 m/s at a = 2 deg yielded 
a change in Cn/Ca of less than 0.00015, so that the 
predicted ratio of normal to axial acceleration has an 
uncertainty of at least ±0.00015. The change in C m 
for this refinement was less than 0.0004. In a related 
grid coarsening test the solution at the 6000 m/s and 
a = 2 deg. was computed on a grid with every other 
mesh point removed in all three computational direc- 
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Fig. 9 Reconstructed density and ratio of mea- 
sured accelerations as function of velocity for Mars 
Pathfinder during flight through the continuum, 
hypersonic flow regime. 

tions. This test yielded a change in Cn/Ca of less 
than 0.00012 and a change in C m less than 0.0003. 

Measured axial accelerations in the hypersonic, con- 
tinuum domain are accurate to within 1%. Uncertain- 
ties in measured normal accelerations are more difficult 
to characterize; they are expected to be less than 1% 
of full scale (A An < .008g) in this domain. Velocities 
are derived by integrating the measured accelerations 
in time from the initial state vector. Uncertainties in 
the derived velocities should be less than 1% because 
random errors in measured accelerations smooth out 
and calibration checks remove bias errors in measured 
accelerations. Furthermore, an independent check of 
derived velocity with altimeter data prior to parachute 
deployment showed differences less than 1%. 

Uncertainties in the derived density are linearly re- 
lated to the total uncertainty in the simulated values 
Ca which include uncertainties associated with ar and 
gas chemistry model. (Recall Eq. 2.) Uncertainties in 
the derived density are quadratically related to the 
total uncertainty in freestream velocity. Comparison 
of LAURA simulations to ground based experimental 
data for an axisymmetric blunt body (METEOR 5 ) at 
Mach 6 in which uncertainties in gas chemistry (perfect 
gas) and ar (sting mounted) are near zero (relative to 
the flight case) show differences less than 2% in Ca, 
and less than 5% in Cjv and C m . Uncertainty in Ca 
associated with the gas model is difficult to define, 
but thought to be approximately 1% for low angles 
of attack based on computational analysis in the next 
section. Uncertainty in the derived ar is implicitly 
related to the uncertainty in Cjs and measured accel- 
erations; and, given the values discussed above and the 
tables and figures from the preflight analysis, 1 the un- 


certainty in ax should be less than 0.6 degrees. The 
uncertainty in Ca due to the uncertainty in cut is less 
than 0.7%. Effects of base pressure on the hypersonic 
aerodynamic drag are less than Q.5%. 1 Consequently, 
the total RMS value for uncertainty in Ca in the hy- 
personic, continuum domain at low angles of attack 
is less than 3%. The corresponding total RAIS value 
for uncertainty in p ^ in this domain is less than 4%. 
Comparison of computed drag coefficient 1 to experi- 
mental data 6 for the Viking configuration in ballistic 
range tests in CCU fall within this 3% estimate. 

As a final note, all of the simulations assume lami- 
nar flow. The assumption should be valid for most of 
the hypersonic, continuum trajectory presented here. 
At the Ex, = 3996 m/s trajectory point ( M ^ = 18, 
Reo , oo = 1-2 million), the value of Reg/M e is less 
than 500 over most of the cone and drops below zero as 
flow expands over the shoulder. In recent comparisons 
to engineering code assessments of Reg/M e the CFD 
valuation of this quantity is roughly twice the engi- 
neering code valuation. Assuming a transition criteria 
of Reg/M e = 250 for blunt bodies from engineering 
code assessments turbulent flow may occur near this 
condition. The contribution of laminar skin friction 
to Ca is less than 0.2% as noted in block 2 of Table 

4 (Moo = 16, Uoo = 3515 m/s, a = 2 deg.) in the 
preflight analysis. 1 If we assume: (1) no relaminariza- 
tion at the shoulder where the strongest contribution 
of shear to Ca will occur; (2) negligible displacement 
thickness effects of a turbulent boundary layer on pres- 
sure for this very blunt configuration; and (3) a factor 

5 increase in shear over laminar values, then an ad- 
ditional uncertainty of 1% in Ca is introduced. The 
contribution from turbulent shear at the shoulder to 
C m at small angles of attack is more difficult to esti- 
mate because of the greater importance of differences 
between windside shoulder and leeside shoulder con- 
tributions. The laminar shear contribution to C m for 
the Moo = 16 simulation noted above is a destabilizing 
increment equal to 0.00012. Again assuming a factor 
5 increase in the turbulent contribution leads to an 
additional uncertainty in C m equal to 0.0006. 

Sensitivity Analysis 

The aerodynamic data base used in the POST code 
was originally generated on the basis of a simulated 
trajectory for a 418 kg vehicle with a ballistic coef- 
ficient of 45 kg/m 2 in the Mars-GRAM atmospheric 
model. 7 The vehicle grew as the design matured and 
the actual entry mass was 585.3 kg with a ballistic co- 
efficient of 62 kg/m 2 . In the final design phase, cruise 
and entry operations, and post-flight analysis, the 
Clancy atmospheric model 8,9 was used. The POST 
code used velocity as a curve fit parameter for aero- 
dynamic coefficients above Mach 12 in the continuum 
regime because of the important role of effective 7 on 
aerodynamics and its dependence on total enthalpy. 
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Below Mach 8, Mach number was used as a correlat- 
ing parameter to emphasize compressibility effects on 
aerodynamics. Between Mach 8 and 12, a bridging 
function is used. As a result of differences in mass, 
entry angle, and atmospheric model, a comparison of 
the early preflight and post-flight trajectories will yield 
different densities for the same velocity. As a point of 
reference, the data base used aerodynamic coefficients 
simulated at Foo = 3996 m/s and p^ = 1.55 10 -3 
kg/m 3 ; the density at that same velocity in the recon- 
structed atmosphere was p^ = 1.13 10 -3 kg/m 3 . 

In the continuum, hypersonic flight regime, the 
change in density would be expected to influence aero- 
dynamic coefficients through Reynolds number and 
gas kinetics. As discussed in the previous section, the 
influence of Reynolds number is expected to be small, 
but the influence of gas kinetics on the effective 7 is 
unknown. The primary purpose of this section is to 
quantify the sensitivity of the aerodynamic data base 
to perturbations in density and the gas kinetic model. 
A secondary purpose of this section is to critically eval- 
uate the transition regions between statically stable 
and unstable regions and check that the simulation is 
consistent with measured normal accelerations. 

The moment coefficient at a = 2 deg. for the 
postflight, simulated trajectory is plotted over the sim- 
ulated Pathfinder measurements of normal to axial 
acceleration as a function of velocity in Fig. 10. The 
value of C m at a = 2 deg., designated C m { 2), is 
proportional to C m>a at a = 0 assuming a forward dif- 
ference approximation C mm { 0) = (C m ( 2) — 0) * 90/7 r. 
The purpose of this plot is to identify values of C m ( 2) 
that are associated with the transition between regions 
of static instability and static stability. The transition 
region is delineated with the rectangles. 

We introduce an assumption that critical values in 
C m { 2 ) may be defined that serve to delineate the var- 
ious stability regions as determined from measured 
accelerations. In each transition region an upper limit 
on C m ( 2), C m x r , is identified such that any simulation 
with C m ( 2) > C m< u would be expected to be in a re- 
gion of static instability. In like manner, a lower limit 
on C m { 2), C m< L, is identified such that any simulation 
with C m (‘2 ) < C m _L would be expected to be in a re- 
gion of static stability. If C m _u > C m (2 ) > C TOj l one 
would expect measured accelerations to be within the 
transition between the two stability regions. 

Values of C m _u and C m L are obtained from analysis 
of four different simulations of the Pathfinder trajec- 
tory through the reconstructed atmosphere. The ve- 
hicle was numerically flown using the 6-DOF POST 
code with the baseline aerodynamic data base and 
three perturbations in which the curve fits of aero- 
dynamic coefficients as a function of velocity were 
shifted by +1%, -5%, and -10% in velocity space (i.e. 
C x = CxiVoo — O.lFoo)). This shift causes the instabili- 
ties to arise earlier or later in the simulated trajectory 



Fig. 10 Ratio of normal to axial accelerations and 
moment coefficient at a = 2 deg. for the preflight 
simulated trajectory as a function of velocity. The 
transition between regions of static instability and 
static stability are bounded with the rectangular 
boxes. 

and allows us to check the assumption that critical 
values for C m 0 and C, n> i can be defined. Values for 
C mt u and C m ,L based on this subjective selection crite- 
rion are respectively -0.0003 and -0.0013 (see bounding 
boxes in Fig. 10). When the major contributor to un- 
certainty in C m , CG location, is included the transition 
range expands to +0.00034 to -0.00194. 


Results and Discussion 

Four test points are chosen along the reconstructed 
trajectory and are identified in Table 1 and in Fig. 
11. The simulated data in Fig. 11 is based on flight 
through the reconstructed atmosphere, in contrast to 
the data in Fig. 7 a) based on flight through the 
preflight atmospheric model. Except for the lack of 
high frequency oscillations in the measured data be- 
tween 3400 m/s and 6000 m/s there is fairly good 
qualitative agreement between the measured and simu- 
lated data. (Temporal resolution of the accelerometers 
should have picked up such oscillations if they were 
present.) The magnitude of measured accelerations is 
slightly larger (15-30%) than the temporally averaged 
level of simulations in the statically stable region be- 
tween 4500 and 6000 m/s and in the statically unstable 
region between 3400 and 4000 m/s. The measured du- 
ration of the statically stable region in velocity space 
is shorter than the simulated duration. These obser- 
vations are consistent with trends shown in Fig. 9 of 
Braun et. al. 2 in which the CG moves off the axis. 
The current simulation assumes a 0.448 mm offset CG 
location. 
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Fig. 11 Comparison of An / Aa as measured in 
flight and as simulated with the POST code in the 
reconstructed atmosphere using the CFD derived 
data base for aerodynamic coefficients. 



Fig. 12 The computed sensitivity of Pathfinder 
aerodynamic coefficients on freestream density for 
a velocity of 6000. km/s. 


Table 1 - Test Points and Baseline Aerodynamics 


Pt. 

V — 

^ OO q 

kg 

Poo? 7n3 

C A ( 2) 

C N ( 2 ) 

C m ( 2) 

A 

6000 

4.53 10 “ 4 

1.707 

0.004399 

-0.001292 

B 

4270 

1.00 10“ 3 

1.714 

0.004394 

-0.001929 

C 

3996 

1.13 10 -3 

1.712 

0.003250 

-0.000812 

D 

3648 

1.26 10 -3 

1.698 

0.001883 

+0.000331 


The aerodynamic coefficients for the baseline eight 
species (COo, CO, No, Oo, NO, C, N, 0) kinetic 
model 10 are included in the Table 1. Point A is in the 
middle of the transition region, near peak heating and 
prior to peak dynamic pressure where the sonic line 
is completing its movement to the nose with decreas- 
ing 7 . Point B is in the middle of the next transition, 
after peak dynamic pressure in which the sonic line 
is beginning to move back to the shoulder from the 
nose with increasing 7 associated with decreasing flow 
enthalpy. At this point in Fig. 11 the POST simula- 
tion indicates a region of static stability, in apparent 
disagreement with the flight data. Point C is at the 
end of this transition region containing Point B; the 
velocity here corresponds to one used in establishing 
the aerodynamic data base. Point D is at the local 
maximum for A n /Aa ■ 

Given the assumptions stated previously, values of 
C m ( 2) for Points A-C would be expected to be be- 
tween C mt u and C m ^L- The value of C m ( 2) for Point 
D would be expected to be greater than C m ^u- As is 
evident in Table 1, the baseline kinetic model yields 
values of C m ( 2 ) which are consistent with the critical 
values of C ra ,r/ and C mt L for points A, C, and D, even 
without adding the uncertainty bar associated with 
CG location. The uncertainty bar around point B as- 
sociated with CG location just overlaps the lower limit 



Fig. 13 The computed sensitivity of Pathfinder 
aerodynamic coefficients on dissociation rate coef- 
ficient for CO 2 at a freestream density 0.0004532 
kg/m 3 and velocity of 6000. km/s. 

for C m< L that defines the transition region. 

The sensitivity of (7 4 (2) and C m { 2) with respect to 
Poo and kinetic model at point A are shown respec- 
tively in Figs. 12 and 13. Open symbols represent 
the baseline grid results (60 x 30 x 64) and the filled 
symbols show a result with a factor two coarsening 
in all coordinate directions. Aerodynamic coefficients 
are insensitive to perturbations as large as 10 % about 
the nominal density of 4.53 10 -4 kg/m 3 . The value 
for C m ( 2) at Point A is at the lower edge of the 
expected transition range. Aerodynamic coefficients 
are also insensitive to changes in the dissociation rate 
constant for CO 2 by factors varying from 0.5 to 2.0. 
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Fig. 14 The computed sensitivity of Pathfinder 
aerodynamic coefficients on freestream density for 
a velocity of 3996. km/s. 

Inspection of CCD mass fraction in all three solutions 
shows that CCD is substantially depleted, other species 
play a more dominant role in establishing 7 . Chang- 
ing the effective temperature for dissociation in all 
such reactions at Point A from the baseline model 
(T eff = T°- 5 Ty 5 , open symbols) to T eff = T 07 Tf 3 
(filled symbols) has almost no effect. While thermal 
nonequilibrium is present in the shock layer there is 
not enough thermal nonequilibrium over substantial 
portions of the shock layer to make a significant effect 
on aerodynamics through effective kinetic rates. 

The sensitivity of 6 ) 4 (2) and C m { 2) with respect 
to Poo and kinetic model at point C are shown respec- 
tively in Figs. 14 and 15. Aerodynamic coefficients are 
still insensitive to perturbations as large as 40% about 
the nominal density of 1.13 10 -3 kg/m 3 . In contrast, 
the effect on C m ( 2) of a 7% perturbation in velocity 
at the same density is large, as observed by comparing 
point B in Table 1 with the result for equal to .001 
kg/m 3 in Fig. 14. The effects of dissociation rate con- 
stant for CO 2 on aerodynamics is stronger at point C 
than at Point A. Previous characterization of the flow 
at this point as “near equilibrium” over-simplified the 
role of gas kinetics during this portion of the trajec- 
tory. Doubling the rate factor moves C m ( 2) closer to 
a predicted domain of static stability. Reducing the 
rate by a factor of 10 moves C m ( 2) closer to a pre- 
dicted domain of static instability. Nevertheless, the 
change in Ca( 2 ) over this entire range of kinetic model 
adjustments is less than 0 . 6 %. 

The movement of the sonic line associated with 
changes in the CO 2 dissociation rate at Point C are 
shown in Fig. 16. The baseline kinetic model results 
are in Fig. 16 c, with rate factor equal to 1. Increas- 
ing rate factor causes the sonic bubble to shrink and 


Fig. 15 The computed sensitivity of Pathfinder 
aerodynamic coefficients on dissociation rate coef- 
ficient for CO 2 at a freestream density 0.001127 
kg/m 3 and velocity of 3996. km/s. 



(a) 0.1 (b) 0.5 (c) 1 .0 (d) 2.0 


Fig. 16 The computed sensitivity of Pathfinder 
sonic line location on dissociation rate coefficient 
for CO 2 at a freestream density 0.001127 kg/m 3 
and velocity of 3996. km/s. 

the moment coefficient to move into the domain of 
static stability. Decreasing rate factor causes the sonic 
bubble to grow and approach a condition of static in- 
stability. The sonic bubble emanates from an inflection 
point in the bow shock which eventually disappears 
with the merging of the bubble with the sonic line on 
the edge of the boundary layer. 

The sensitivity tests at Points A and C show that 
aerodynamic database is insensitive to perturbations 
in density as large as 10%. It is believed this lack of 
sensitivity would persist across the entire hypersonic, 
continuum domain for the Pathfinder vehicle shape. 
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The sensitivity will manifest itself during the transi- 
tional flow regime prior to continuum; in fact, Moss 
et. al. 11 predicted static instabilities in the transitional 
regime for Pathfinder as a function of Knudsen number 
variation. 

The baseline thermochemical model yields values for 
C m (2) which are consistent, within estimated uncer- 
tainty bands, of the measured normal accelerations. 
Point B shows the greatest disagreement between 
measurement and simulation of normal accelerations. 
Agreement would likely be improved using a simula- 
tion with a larger offset CG and/or slower CO 2 disso- 
ciation rates to define C m . 

The normal acceleration data is useful as a con- 
sistency check for the simulation models, much like 
measurements of shock standoff distance in ground 
based experiments are used. While one could not use 
the data by itself to propose changes to any element (s) 
of the chemical kinetic model, one should question 
any model changes that cause inconsistent values of 
C m {2 ) with measured normal accelerations along any 
point on the Pathfinder trajectory. In the tests con- 
ducted herein, rather large perturbations in the kinetic 
model are required before computed values for C m ( 2) 
are pushed outside of the expected limits based on 
measured normal accelerations. 

Conclusions 

(1) Predictions of two regions of static instability 
during the hypersonic, continuum portion of the Mars 
Pathfinder trajectory are verified by flight measure- 
ments of normal accelerations. The relation of these 
instabilities to sonic line movement between the nose 
and the shoulder as a function of trajectory point and 
kinetic model is demonstrated. 

(2) Perturbations to density as large as 10% about 
nominal values for the reconstructed trajectory dur- 
ing the hypersonic, continuum regime cause negligible 
change (less than 1%) in Ca ■ Large changes in the 
kinetic model in this regime cause changes in Ca of 
approximately 1%. Given other sources of uncertainty 
unrelated to the simulation model, the reconstructed 
density, with uncertainty linearly related to uncer- 
tainty in Ca, is expected to have approximately 4% 
uncertainty in the hypersonic, continuum regime. 

(3) The occurrence of static instabilities as recorded 
by normal accelerations during the Mars Pathfinder 
entry, descent, and landing (EDL) provide an opportu- 
nity for computational fluid dynamic (CFD) code val- 
idation. Computed moment coefficients are consistent 
with identified ranges for measured normal accelera- 
tions indicating flight in a statically stable or unstable 
regime. The definitions of these ranges for C m are 
independent of CFD simulation; they are extracted 
from 6 degree-of-freedom trajectory simulations us- 
ing several aerodynamic coefficient models for flight 
through the Mars atmosphere. The present simula- 


tions are consistent with measured data within the 
defined ranges, including the uncertainty band. 
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